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We present here a intra-nuclear cascade model implemented in Geant4 5.0. The cascade model is based on 
re-engineering of INUCL code. Models included are Bertini intra-nuclear cascade model with exitons, pre- 
cquilibrium model, nucleus explosion model, fission model, and evaporation model. Intermediate energy nuclear 
reactions from 100 MeV to 3 GeV energy are treated for proton, neutron, pions, photon and nuclear isotopes. 
We represent overview of the models, review results achieved from simulations and make comparisons with 
experimental data. 



1. INTRODUCTION 

The intra-nuclear cascade model (INC) was first 
proposed by Serber in 1947 [lj- He noticed that, in 
particle-nuclear collisions the deBroglie wave-length of 
the incident particle is comparable to or shorter than 
the average intra-nucleon distance. Hence, the jus- 
tification for describing the interactions in terms of 
particle-particle collisions. 

The INC has been successfully used in the Monte 
Carlo simulations at intermediate energy region since 
Goldberger made first calculations by hand in 1947 @ ■ 
First computer simulations were done by Metropolis 
et al. in 1958 [§}. Standard methods in INC imple- 
mentations were formed when Bertini published his 
results in 1968 J4J. An important addition was exci- 
ton model introduced by Griffin in 1966 [T^ . 

Our presentation describes the implementation of 
a Bertini INC model in the Geant4 hadronic physics 
frameworks 0. Geant4|2| is a Monte Carlo particle 
detector simulation toolkit, having applications also 
in the medical field and space sciences. Geant4 pro- 
vides a flexible framework for modular implementa- 
tion of various kinds of hadronic interactions. Geant4 
exploits advanced Software Engineering techniques 
and Object Oriented technology to achieve the trans- 
parency of the physics implementation and to this way 
provide the possibility of validating the physics results 
more easily. 

At lcvel-2, the hadronic model frameworks are 
based on concepts of physics processes, cross-section 
and models for final state generation. While the pro- 
cess is a general concept, models are allowed to have 
restrictions in process type, material, element and en- 
ergy range. Several models can be utilized by one 
process class; for instance, a process class for inelastic 
scattering can use distinct models for different ener- 
gies. 

Process classes utilize model classes to determine 
the secondaries produced in the interaction and to cal- 
culate the momenta of the particles. Here we preset 
a collection of such models providing medium-energy 
intra-nuclear cascade treatment. 



2. GEANT4 CASCADE MODEL 

In inelastic particle-nucleus collisions a fast phase 
(10~ 23 - 10~ 22 s) of INC results in a highly exited nu- 
cleus, and is followed by fission and pre-equilibrium 
emission. A slower (1CP 18 — 10~ 16 s) compound nu- 
cleus phase follows with evaporation. A Boltzman 
equation must be solved to treat the physical process 
of the collision in detail. 

The intra-nuclear cascade model developed by 
Bertini Q- 0- solves the Boltzmann equation on the 
average. The model has been implemented in several 
codes such as HETC Q. Our model is based on re- 
engineering of INUCL code ||. Models included are 
a Bertini intra-nuclear cascade model with exitons, 
a pre-equilibrium model, a simple nuclear explosion 
model, a fission model, and a evaporation model. 

The nuclear model consist of a three-region approx- 
imation to the continuously changing density distri- 
bution of nuclear matter within nuclei. Relativistic 
kinematics is applied throughout the cascade. The 
cascade is stopped when all the particles which can 
escape the nucleus, have done so. Then conformity 
with the energy conservation law is checked. 

2.1. Model limits 

Particles treated are nucleons, pions, photons, and 
nuclear isotopes. A bullet particle can be a proton, 
neutron or pion. The range of targets allowed is arbi- 
trary. 

The necessary condition of validity of the INC 
model is Xb/v « t c « At, where \b is the de- 
Broglie wavelength of the nucleons, v is the aver- 
age relative nuclcon-nucleon velocity and At is the 
time interval between collisions. The physical founda- 
tion becomes approximate at energies less than about 
200MeV, and there needs to be supplemented with a 
pre-equilibrium model. Also, at energies higher than 
5-10 GeV the INC picture breaks down. Our imple- 
mentation in Geant4 has been tested with bullet ki- 
netic energy between 100 MeV and 3 GeV. 
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400 MeV 




Figure 1: Schematic presentation of the intra- nuclear 
cascade. A hadron with 400 MeV energy is forming an 
INC history. Crosses present the Pauli exclusion 
principle in action. (The picture is a reproduction from 
original work of Bertini Q.) 

2.2. Intra-nuclear cascade model 

The basic steps of the INC model are summarized 
below: 

1. The spatial point, where the incident particle 
enters, is selected uniformly over the projected 
area of the nucleus. 

2. Total, free particle-particle cross-sections and 
region-dependent nucleon densities are used to 
select the path length for the projectile particle. 

3. The momentum of a struck nucleon, the type of 
reaction, and the four momentum of the reaction 
products are determined. 

4. the exciton model is updated as the cascade pro- 
ceeds. 

5. If Pauli's exclusion principle allows and 
Epartide > E cuto ff = 2 MeV, step (2) is per- 
formed to transport the products. 

After INC, the residual excitation energy of the re- 
sulting nucleus is used as input for a non-equilibrium 
model. The schematic presentation of Bertini's INC 
is shown in Fig. ^ 

2.3. The nuclear model 

Some of the basic features of the nuclear model are 
the following: 

• The nucleons are assumed to to have a Fermi 
gas momentum distribution. The Fermi energy 
is calculated in a local density approximation i.e. 
it is made radius dependent with Fermi momen- 
turn p F (r) = (^Mli)i 



• The nucleon's binding energies (BE) are calcu- 
lated using a mass formula. A parametrization 
of the nuclear binding energy uses a combination 
of the Kummel mass formula, and experimental 
data. Also, a asymptotic high temperature mass 
formula is used, for cases where it is not possible 
to use experimental data. 

2.3.1 . Initialization 

The initialization phase fixes the nucleus radius and 
momentum according to the Fermi gas model. 

If the target is Hydrogen (A = 1), a direct particle- 
particle collision is performed, and no nuclear model- 
ing is used. 

If 1 < A < 4, a nuclei model consisting one layer 
with radius of 8.0 fm is created. 

For 4 < A < 11, a nuclei model is composed of three 
concentric spheres i = {1,2,3} with radii 



n( ai ) = ^C?(l - -) + 6Ay/-log( ai ) 

where a* = {0.01,0.3,0.7} and Ci = 3.3836.4 1 / 3 
If A > 11, nuclei are modeled with three concentric 
spheres as well. The sphere radii are then defined as: 

r<(aO = C 2 log( 1 + e - 1) + d 

Cti 

where C 2 = 1.7234. 

The potential energy for nucleon N is 

2 

V N = ^ + BE N (A,Z) 

where pp is a Fermi momentum and BE a binding 
energy. 

Impulse distribution in each region follows Fermi 
distribution with zero temperature. 

f(p) = cp 2 (1) 

where 

i-pf 

/ f(p)dp = n p or n n . (2) 
Jo 

Here n p and n n are the number of protons or neu- 
trons in a region and P f is momentum corresponding 
the Fermi energy 

^ = ^ = &i— ) § > (3) 
Ittin IrriN v 

which depends on the density n/v of particles, and 
which is different for each particle and each region. 
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2.3.2. Pauli's exclusion principle 

Pauli's exclusion principle forbids interactions 
where the products would be in occupied states. Fol- 
lowing the assumption of a completely degenerate 
Fermi gas, the levels are filled from the lowest level. 
The minimum energy allowed for a collision prod- 
uct corresponds to the lowest unfilled level of sys- 
tem, which is the Fermi energy in the region. So in 
practice, Pauli exclusion principle is taken into ac- 
count by accepting only secondary nucleons which 
have E N > Ep. 

2.3.3. Cross-sections and kinematics 

Path lengths of nucleons in the nucleus are sam- 
pled according to the local density and free nucleon- 
nucleon cross-sections. Angles after collisions are 
sampled from experimental differential cross-sections. 
Tabulated total reaction cross-sections are calculated 
by Letaw's formulation [Ti El Il6|. For nucleon- 
nucleon cross-sections, parameterizations based on the 
experimental energy and isospin dependent data. The 
parameterization described in ^3 is used. 

For pion the INC cross-sections are provided to 
treat elastic collisions, and inelastic channels: 7r~n — > 
7T°n, 7r°p — > 7r + n and 7r°n — > 7r~p. Multiple particle 
production is also implemented. 

The S-wave pion absorption channels 7r + nn — > pn, 
7r + pn — > pp, 7r°nn — > X , 7r n pn — > pn, 7r°pp — ► pp, 
7r~nn — > X , 7r _ pn — > nn , and 7r~pp — > pn are imple- 
mented. 

2.4. Pre-equilibrium model 

The Geant4 cascade model implements the exciton 
model proposed by Griffin [l^. In this model nu- 
clear states are characterized by the number of exited 
particles and holes (the exitons). INC collisions give 
rise to a sequence of states characterized by increasing 
exciton number, eventually leading to a equilibrated 
nucleus. For practical implementation of the exciton 
model we use level density parameters from .18) and 
the matrix elements from [l9j . 

In the exciton model the possible selection rules for 
particle-hole configurations in the course of the cas- 
cade are: Ap = 0, ±1 Ah = 0, ±1 An = 0, ±2, where 
p is the number of particle, h is number of holes and 
n = p + h is the number of exitons. 

The cascade pre-equilibrium model uses target ex- 
citation data, and exciton configurations for neutron 
and proton to produce the non-equilibrium evapora- 
tion. The angular distribution is isotropic in the frame 
of rest of the exciton system. 

The parameterizations of the level density used, are 
tabulated both with their A and Z dependence and 
including a high temperature behavior. The nuclear 
binding energy is using a smooth liquid high energy 
formula. 



2.5. Break-up models 

Fermi break-up is allowed only in some extreme 
cases, i.e. for light nuclei (A < 12 and 3(A — Z) < Z < 
6 ) and if E excit ation > 3E binding ■ A simple explosion 
model decays the nucleus into neutrons and protons 
and decreases exotic evaporation processes. 

The fission model is a phenomenological model 
using potential minimization. Binding energy 
parametrization is used and some features of the fis- 
sion statistical model are incorporated as in[20|. 

2.6. Evaporation model 

The statistical theory for particle emission from ex- 
ited nuclei remaining after INC was originally devel- 
oped by Weisskopf $] ■ This model assumes complete 
energy equilibration before particle emission, and re- 
equilibration of excitation energies between successive 
evaporation emissions. As a result, the angular distri- 
bution of emitted particles is isotropic. 

The Geant4 evaporation model for cascade im- 
plementation adapts the widely used computational 
method developed by Dostrowski [IcL lllj . The emis- 
sion of particles is computed until the excitation en- 
ergy fall below a cutoff. If a light nucleus is highly 
exited, a Fermi break-up model is executed. In addi- 
tion, fission is performed when the fission channel is 
open. The main chain of evaporation is followed until 
Excitation falls below E cuto ff = 0.1 MeV. The evap- 
oration model ends with a 7 emission chain, which is 
followed until E excitaUo n < Kutoff = 10 ~ 15 MeV - 

3. IMPLEMENTATION 

The Bertini cascade model is implemented in 
the Geant4 hadronic physics framework. The 
Bertini cascade source code module is located in di- 
rectory Geant4/ source/processes /hadronic/models / - 
cas cade /cascade. All the models are used collectively 
through interface method Apply Yourself defined in a 
class GJ^Cascadelnterface. A Geant4 track (G4Track) 
and a nucleus (G4Nucleus) are given as a parameters. 

Currently implementation is not optimized for 
speed. The typical timing behavior for different target 
nuclons and bullet energies is presented in Fig. [5] 



4. RESULTS 

We have tested the physics performance of Bertini 
cascade models in the first Geant4 5.0 implementa- 
tion for proton induced reactions for energies from 
100 MeV to 3 GeV. Detailed comparisons with ex- 
perimental data has been made in energy range 160 - 
800 MeV. 
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Figure 2: Speed of single INC performed for Al, Fe and 
Pb nuclei for proton bullet at energy-range 1-8 GeV. A 
single PC with 300 MHz PII is used. 
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Figure 4: p(256MeV) + Fe -> n(0 = 30°) + X. 



Figure 3: Geant4 Bertini INC model simulation of 
p(256MeV) + Fe -> n(9 = 7.5°) + X in comparison with 
experimental data [2l| . 
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The double differential cross-sections for neutrons 
at angles of 7.5°, 30°, 60°, and 150° are presented in 
Figs. 13161 Proton incident particle at 256 MeV energy 
is hitting an iron target and all Bertini sub-models are 
used. Data were taken from ,21]. 

Figures [7| and [S] show double differential cross- 
section of neutrons resulting from 597 MeV protons 
impinging on lead. Data were taken from 

Finally, an example of pion production physics per- 
formance. Double differential cross-sections for tt + are 
given for angels from 22.5° to 135° in Figs.EH3 

The agreement with experimental data is found to 
be from moderate to relatively good. Data fo the 
585 MeV pion comparisons were taken from |2^| . 



Figure 5: p(256MeV) + Fe -» n(9 = 60°) + X. 

5. CONCLUSION 

We have released a Bertini INC model in Geant4. 
Exitons, pre-equilibrium, nucleus explosion, fission, 
and evaporation are modeled. Particles treated are 
7, 7T, n, p, and nuclear isotopes, and incident nucle- 
ons and pions can be treated. We have tested the 
code in a energy range 60 MeV - 3 GeV and found a 
reasonable agreement with experimental data. Tun- 
ing of performance and models are planned for future 
Geant4 releases. We also plan to extend the upper 
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Figure 6: p(256MeV) + Fe -> n{6 = 150°) + X. 
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Figure 7: Geant4 Bertini INC model simulation of 
p(597MeV) + Pb -> n{6 = 60°) + X in comparison with 
experimental data [2^|. 



energy limit of the models. 



References 

[1] R. Serber, "Nuclear Reactions at High Energies" , 

Phys. Rev. 72, 1948 
[2] M. Goldberger, "The Interaction of High Energy 

Neutrons and Heavy Nuclei" , Phys. Rev. 74, 1948 




i I i I i I I I U U I 

100 200 300 400 

EuJMeVJ 
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Figure 11: p(585MeV) + Pb -> n + (d = 60°) + X. 
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